Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  546 1691 (N20E53 L= 22)   706    187        2 ( 518/1647- 578/1741) -33.7 E58-E48 N17-N23
      2  308 1588 (N09E72 L=  3)   863    123        3 ( 279/1556- 341/1622)  68.2 E74-E68 N06-N12
      3  418 1527 (N07E55 L= 20)  1934    406        7 ( 340/1440- 520/1617)  58.4 E59-E50 N00-N14
      4 1297 1349 (N19W06 L= 81)   888    257        1 (1271/1229-1329/1469) -68.2 W03-W10 N12-N27
      5 1630 1344 (N30W29 L=104)   298     21        0 (1609/1325-1653/1362)   0.0 W26-W30 N29-N31
      6  233 1310 (S09E63 L= 12)  1088    184        2 ( 177/1276- 271/1351)  35.5 E69-E58 S11-S06
      7 1486 1181 (N15W22 L= 97)   404    114        0 (1467/1119-1504/1240) -60.3 W19-W23 N11-N18
      8  192  950 (S31E57 L= 18)  1680    429        2 ( 101/ 923- 284/ 982)  14.7 E66-E47 S33-S27
      9 1306  909 (S05W19 L= 94)   230     42        2 (1284/ 879-1327/ 932) -45.0 W16-W19 S06-S03
     10  750  874 (S22E12 L= 63)   253     76        1 ( 728/ 831- 775/ 912)  90.0 E12-E10 S25-S19
     11 1853  857 (N10W61 L=136)   548     85        2 (1831/ 825-1870/ 895) -53.1 W58-W62 N08-N12
     12  464  739 (S38E29 L= 46)  1371    325        3 ( 352/ 707- 588/ 791)  16.7 E39-E19 S42-S33
     13 1837  667 (N00W69 L=144)  1330    239        1 (1803/ 642-1888/ 707)  10.6 W64-W80 S01-N03
     14  966  590 (S32W10 L= 85)   884    201        1 ( 883/ 522-1052/ 676)  70.3 W06-W12 S38-S24
     15 1307  579 (S22W31 L=106)   579    165        1 (1223/ 543-1371/ 639)  66.0 W26-W33 S26-S17
     16 1571  579 (S14W48 L=123)   641    209        2 (1506/ 547-1631/ 624)   0.0 W41-W54 S15-S11
     17 1317  436 (S30W39 L=114)   530    109        0 (1272/ 413-1369/ 479)  63.4 W36-W40 S31-S25
  

    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