Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1509 1565 (N30W39 L=132)  1823    365        4 (1356/1456-1641/1671) -41.4 W30-W47 N22-N38
      2 1696 1590 (N31W61 L=154)  1956    170        1 (1634/1523-1763/1634) -45.0 W54-W67 N26-N34
      3  787 1478 (N27E14 L= 79)   267     80        0 ( 762/1437- 816/1508)  51.3 E15-E11 N24-N29
      4 1077 1433 (N23W06 L= 99)   231     91        0 (1050/1396-1102/1480) -71.6 W04-W06 N20-N26
      5  359 1369 (N23E46 L= 47)   882    160        5 ( 294/1321- 415/1425)  33.7 E51-E41 N20-N26
      6 1568 1382 (N17W39 L=132)   593    102        0 (1505/1353-1616/1407) -14.0 W34-W42 N15-N19
      7 1236 1357 (N17W15 L=108)   455    160        1 (1173/1332-1305/1386)   0.0 W11-W19 N15-N18
      8 1232  925 (S09W12 L=105)   603    196        1 (1166/ 864-1307/ 996)  33.7 W07-W17 S12-S05
      9 1199  688 (S24W10 L=103)   487    138        3 (1178/ 616-1217/ 764)  90.0 W08-W10 S28-S18
     10 1715  654 (S28W51 L=144)   575     63        1 (1682/ 624-1738/ 681)  26.6 W47-W54 S29-S25
     11 1635  484 (S39W50 L=143)   952    105        3 (1610/ 428-1662/ 542) -63.4 W45-W54 S43-S35
     12  702  299 (S49E37 L= 56)   433     37        1 ( 682/ 274- 720/ 326)  53.1 E40-E34 S50-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