Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  604 1797 (N32E61 L=264)  1044    107        1 ( 566/1752- 641/1835) -12.8 E74-E52 N29-N34
      2  518 1629 (N18E52 L=273)   895     87        4 ( 488/1583- 551/1678) -50.2 E55-E49 N15-N21
      3 1332 1444 (N28W07 L=332)  1285    443        0 (1254/1262-1410/1616) -46.6 E02-W16 N18-N37
      4 1562 1350 (N29W25 L=350)   644    118        0 (1512/1300-1608/1384) -14.0 W21-W29 N27-N30
      5  959 1202 (N04E08 L=317)   731    225        4 ( 857/1149-1015/1281)  53.1 E12-E05 N00-N09
      6  259 1232 (S10E57 L=268)   495     68        2 ( 237/1200- 294/1257)  63.4 E59-E53 S12-S08
      7 1344 1179 (N12W14 L=339)   242    125        0 (1299/1121-1369/1206) -31.0 W11-W16 N09-N13
      8 1632 1109 (N17W35 L=  0)   442     98        0 (1610/1064-1676/1148) -26.6 W32-W38 N15-N18
      9 1666 1000 (N12W40 L=  5)   283     61        2 (1636/ 976-1689/1030)  63.4 W37-W41 N09-N13
     10 1509  902 (N02W31 L=356)   812    183        1 (1427/ 879-1597/ 940)  14.0 W25-W37 N00-N04
     11 1229  815 (S11W17 L=342)   216     76        0 (1197/ 794-1260/ 853)  63.4 W15-W18 S12-S07
     12  468  806 (S31E31 L=295)   330    125        0 ( 429/ 769- 526/ 830)  40.6 E33-E26 S34-S28
     13 1171  749 (S16W16 L=341)   506    158        1 (1111/ 716-1255/ 778)  35.5 W13-W20 S19-S13
     14 1575  626 (S11W46 L= 11)   332     96        0 (1547/ 589-1598/ 680) -71.6 W44-W46 S13-S07
     15  241  618 (S50E62 L=263)  1244    112        3 ( 182/ 600- 289/ 642)   2.0 E81-E52 S51-S46
     16 1239  521 (S27W29 L=354)   274    121        0 (1189/ 474-1285/ 560) -11.3 W23-W33 S28-S25
     17  968  408 (S42W17 L=342)   543     87        1 ( 927/ 373-1011/ 437)  15.9 W13-W20 S44-S40
     18  690  364 (S56E03 L=322)   305     48        1 ( 670/ 339- 706/ 387) -21.8 E05-W00 S57-S53
  

    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