Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  500 1624 (N16E53 L=189)   380     99        0 ( 475/1589- 516/1650) -33.7 E56-E50 N13-N17
      2  453 1521 (N08E51 L=190)   870    189        0 ( 377/1483- 527/1570)  35.0 E57-E47 N05-N12
      3 1212 1525 (N29E04 L=238)   281     33        0 (1184/1507-1239/1540)  18.4 E05-E01 N27-N30
      4 1304 1386 (N23W06 L=248)   411     66        0 (1284/1350-1326/1430) -68.2 W04-W07 N20-N25
      5  265 1179 (S15E54 L=187)  1364    304        0 ( 173/1095- 326/1304)  61.9 E62-E51 S21-S06
      6 1562 1213 (N21W27 L=268)   250     87        0 (1522/1175-1589/1237)  -9.5 W23-W29 N19-N21
      7 1224 1162 (N07W07 L=249)   279     81        0 (1182/1143-1269/1178)  18.4 W04-W10 N06-N08
      8  664  881 (S23E17 L=225)   449    102        2 ( 639/ 824- 702/ 924)  90.0 E18-E15 S26-S19
      9  182  793 (S41E58 L=184)   583    124        1 ( 153/ 776- 212/ 811)   9.5 E60-E54 S41-S39
     10 1407  710 (S11W31 L=273)   369     60        0 (1368/ 678-1449/ 744)  59.0 W29-W33 S13-S08
     11  970  690 (S26W07 L=249)   290     85        1 ( 935/ 670-1023/ 720)   8.1 W03-W10 S26-S24
     12  785  676 (S32E04 L=238)   327    108        0 ( 718/ 651- 828/ 702)   6.3 E09-W00 S32-S30
     13 1683  505 (S13W62 L=304)   590    129        0 (1659/ 459-1699/ 562) -60.3 W59-W64 S16-S09
     14  811  403 (S49W08 L=250)   599    121        0 ( 791/ 354- 834/ 458) -36.9 W04-W12 S52-S45
     15  945  319 (S50W23 L=265)  1232    194        1 ( 912/ 271-1010/ 395) -26.6 W16-W30 S54-S44
     16 1220  351 (S38W38 L=280)   562    111        2 (1177/ 304-1254/ 391) -63.4 W35-W40 S41-S35
     17 1070  311 (S46W32 L=274)   387     63        0 (1055/ 284-1089/ 343) -45.0 W29-W34 S47-S43
     18  876  287 (S56W21 L=263)   319     58        1 ( 853/ 264- 894/ 309) -11.3 W18-W23 S57-S52
  

    OPENR, u, binfile, /COMP, /GET_LUN
    head = FLTARR(6) ;[xsize,ysize,radius,hh,mm,ss]
    READU, u, head
    data = INTARR(FIX(head[0]),FIX(head[1]))
    READU, u, data
    FREE_LUN, u
    

• each filament has the value of Fil_ID
• all values are in pixel

  c = WHERE(data EQ Fil_ID, area)
  f = '(2(A,I6))'
  PRINT, FORM=f, "Area of Fil_ID", Fil_ID," =", area 

  data = ROTATE(data,7)
  Solar_y = 2048-Solar_y