Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1160 1683 (N39W32 L=143)   439     72        0 (1122/1661-1188/1707)   0.0 W28-W34 N37-N40
      2 1274 1641 (N33W37 L=148)   616    119        0 (1225/1592-1311/1671) -71.6 W33-W38 N28-N35
      3  810 1624 (N46W04 L=115)   286     62        0 ( 776/1611- 846/1638) -18.4 W01-W07 N45-N47
      4 1472 1556 (N22W45 L=156)   437    130        1 (1447/1502-1516/1596)  15.3 W40-W51 N19-N23
      5 1327 1563 (N27W36 L=147)   311     71        1 (1294/1541-1348/1593)   9.5 W32-W38 N24-N28
      6  426 1282 (N34E34 L= 77)   482    142        1 ( 381/1224- 480/1339) -48.4 E38-E29 N29-N38
      7  771 1272 (N25E09 L=102)   326    109        1 ( 748/1225- 800/1323) -81.9 E10-E07 N21-N28
      8  416 1154 (N26E36 L= 75)   286    101        3 ( 379/1124- 462/1181)   0.0 E39-E31 N24-N28
      9  239 1069 (N25E54 L= 57)   679    133        1 ( 203/1044- 292/1111) -11.3 E58-E48 N23-N27
     10 1913 1014 (S21W66 L=177)   631     93        0 (1900/ 968-1925/1058) -80.5 W63-W68 S24-S18
     11  135 1000 (N23E70 L= 41)   962     98        3 ( 106/ 961- 164/1044) -45.0 E75-E64 N20-N25
     12  935  771 (S07E11 L=100)   611    174        1 ( 896/ 684- 986/ 840)  41.2 E15-E06 S11-S04
     13 1709  630 (S39W39 L=150)  1460    314        5 (1575/ 565-1829/ 718)   0.0 W25-W53 S40-S36
     14 1097  506 (S27E10 L=101)   284     65        0 (1063/ 482-1127/ 527)  11.3 E12-E07 S28-S26
     15  976  459 (S27E19 L= 92)   430    113        1 ( 926/ 439-1033/ 483) -33.7 E21-E15 S29-S24
     16  760  404 (S25E35 L= 76)   725    180        1 ( 671/ 379- 845/ 430) -13.0 E41-E28 S26-S23
     17  828  331 (S32E35 L= 76)   383    123        0 ( 785/ 301- 868/ 382)  16.7 E39-E29 S33-S29
     18 1129  191 (S54E26 L= 86)   392     26        0 (1111/ 169-1147/ 209)  23.2 E29-E22 S55-S51
     19 1039  124 (S58E44 L= 67)   638     52        1 (1016/ 104-1057/ 148)  14.0 E50-E37 S60-S54
  

    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