Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1012 1805 (N47E15 L=356)   477     99        1 ( 981/1771-1052/1823)  63.4 E18-E12 N43-N49
      2 1184 1790 (N49E00 L= 11)   432     72        0 (1166/1762-1201/1820) -78.7 E02-W01 N45-N51
      3  892 1704 (N36E20 L=351)   302     87        0 ( 875/1669- 911/1739) -51.3 E22-E18 N33-N38
      4 1476 1599 (N38W25 L= 36)   487     87        1 (1448/1564-1509/1632)  45.0 W22-W27 N34-N40
      5  474 1430 (N13E41 L=330)   258     69        0 ( 439/1407- 493/1449) -21.8 E44-E39 N11-N13
      6  356 1407 (N10E51 L=321)   425    113        1 ( 319/1390- 404/1426)   8.1 E53-E46 N09-N11
      7  403 1298 (N04E44 L=328)   720    172        3 ( 376/1224- 432/1362)  90.0 E46-E41 N00-N07
      8 1891 1196 (N19W61 L= 73)   748    145        1 (1865/1128-1917/1246) -45.0 W57-W64 N15-N21
      9 1854 1169 (N16W57 L= 68)  1099    173        3 (1826/1093-1897/1245) -45.0 W53-W62 N12-N20
     10  327 1133 (S07E48 L=324)   400    111        2 ( 282/1107- 353/1164) -18.4 E51-E45 S07-S04
     11  701 1111 (S05E20 L=351)   238     77        0 ( 663/1075- 731/1144) -31.0 E23-E18 S06-S02
     12 1312  907 (S09W19 L= 30)   470    159        2 (1268/ 838-1352/1001)  79.7 W16-W20 S13-S02
     13  677  948 (S14E20 L=351)   565    159        1 ( 627/ 901- 770/ 989)  33.7 E23-E14 S17-S11
     14  130  898 (S21E67 L=305)  1479    217        2 (  99/ 853- 177/ 954)  33.7 E70-E60 S23-S17
     15 1633  815 (S09W42 L= 54)   251     33        0 (1608/ 801-1655/ 831)  18.4 W40-W43 S10-S07
     16  457  703 (S32E36 L=336)   512     85        2 ( 434/ 660- 485/ 744) -78.7 E37-E33 S34-S29
     17 1677  504 (S26W57 L= 69)   404     57        0 (1650/ 479-1703/ 534)  63.4 W55-W58 S27-S23
  

    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