Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  978 1617 (N33E06 L=293)   637    117        1 ( 948/1561-1015/1661)  81.9 E07-E03 N29-N36
      2  564 1626 (N33E37 L=262)   369     84        0 ( 529/1609- 604/1643)   9.5 E40-E33 N31-N34
      3 1405 1486 (N26W24 L=323)   772    229        0 (1325/1445-1507/1547) -21.0 W18-W31 N23-N29
      4  834 1425 (N19E13 L=285)   198     28        0 ( 817/1398- 848/1448) -76.0 E14-E12 N17-N21
      5 1583 1311 (N16W35 L=334)   829    293        0 (1461/1234-1698/1371) -19.4 W26-W44 N11-N18
      6  221 1329 (N13E60 L=239)   588    101        0 ( 176/1300- 260/1366)  20.6 E65-E56 N11-N15
      7 1885 1193 (N11W63 L=  2)  1665    195        1 (1857/1127-1936/1285) -56.3 W60-W70 N07-N16
      8  170 1246 (N08E64 L=234)   551     64        0 ( 136/1229- 208/1265) -10.3 E70-E59 N07-N09
      9 1600 1039 (S01W36 L=335)   549    150        0 (1575/ 969-1636/1119) -63.4 W34-W39 S04-N04
     10 1173  883 (S12W09 L=308)   201     88        0 (1154/ 845-1205/ 931)  63.4 W08-W11 S14-S08
     11 1841  820 (S11W60 L=359)  1566    131        2 (1791/ 753-1887/ 875)  82.9 W55-W65 S15-S07
     12 1000  676 (S25E00 L=298)  1293    378        0 ( 917/ 559-1165/ 754)  40.2 E05-W09 S32-S19
     13 1656  590 (S27W49 L=348)  1119    158        1 (1573/ 562-1721/ 620)  13.0 W41-W54 S28-S24
     14  262  573 (S33E64 L=235)   723     96        0 ( 226/ 539- 306/ 596)   3.8 E72-E57 S34-S31
     15 1033  536 (S34W03 L=301)   373     70        0 ( 989/ 517-1077/ 557) -21.8 W00-W05 S35-S32
     16 1476  372 (S43W44 L=342)  2676    399        1 (1294/ 346-1658/ 413)   9.7 W25-W64 S45-S38
     17  667  340 (S50E31 L=268)  1838    291        2 ( 507/ 310- 799/ 386)   3.4 E51-E16 S52-S46
  

    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