Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  461 1583 (N13E54 L=162)   680    170        1 ( 440/1519- 480/1651) -56.3 E58-E51 N08-N17
      2 1106 1519 (N25E10 L=206)   239    114        0 (1082/1473-1136/1569)  90.0 E10-E08 N21-N29
      3  741 1373 (N07E26 L=190)   242     48        2 ( 712/1352- 772/1393)  45.0 E27-E24 N05-N08
      4  425 1297 (S05E44 L=172)   612    162        1 ( 401/1206- 468/1349) -83.7 E45-E41 S10-S01
      5  568 1145 (S10E30 L=186)  1676    519        1 ( 397/ 992- 647/1341)  84.8 E39-E24 S22-N00
      6 1374 1232 (N16W14 L=230)   417    166        1 (1318/1166-1401/1306)  84.8 W12-W15 N10-N21
      7 1662 1189 (N23W35 L=251)   672     86        1 (1629/1148-1685/1235) -68.2 W32-W37 N19-N25
      8  136 1100 (S22E67 L=149)  1186    244        1 ( 100/1043- 169/1147)  49.4 E69-E63 S25-S18
      9  369 1094 (S18E43 L=173)   357    126        0 ( 324/1056- 414/1134)  50.2 E46-E40 S20-S14
     10 1267 1000 (S00W14 L=230)   201     58        1 (1246/ 971-1289/1033) -33.7 W12-W15 S01-N01
     11  899  941 (S13E05 L=211)   323    145        0 ( 833/ 899- 978/ 965)  40.6 E07-W00 S17-S10
     12 1787  919 (N12W52 L=268)   322     43        1 (1771/ 891-1803/ 945) -45.0 W50-W53 N10-N13
     13  559  845 (S28E24 L=192)   239     34        1 ( 543/ 823- 581/ 865) -18.4 E25-E22 S29-S26
     14  319  736 (S41E43 L=173)   378     38        1 ( 298/ 711- 342/ 763) -36.9 E45-E40 S42-S39
     15  990  676 (S26W08 L=224)   457    146        1 ( 974/ 609-1020/ 749) -54.5 W05-W11 S29-S21
     16 1161  605 (S25W21 L=237)   405    146        0 (1101/ 586-1234/ 632)  18.4 W16-W25 S26-S21
     17 1093  508 (S32W21 L=237)   521    119        1 (1041/ 458-1128/ 550) -12.5 W16-W25 S34-S30
     18 1658  523 (S13W58 L=274)   421     61        1 (1634/ 500-1680/ 548)  71.6 W55-W61 S14-S11
     19 1147  227 (S48W44 L=260)  1309    222        1 (1114/ 167-1178/ 283) -17.4 W36-W52 S51-S44
  

    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