Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1594 1522 (N16W52 L=339)   896    123        0 (1531/1500-1662/1543) -14.0 W46-W58 N12-N18
      2 1203 1519 (N27W26 L=313)   421     64        0 (1174/1490-1242/1539)   0.0 W23-W29 N25-N28
      3  239 1330 (N43E55 L=231)   702     94        3 ( 199/1300- 289/1355)   0.0 E61-E48 N41-N44
      4 1277 1269 (N12W21 L=307)  1258    455        0 (1112/1161-1444/1329) -43.2 W11-W27 N01-N18
      5 1792 1242 (S05W55 L=342)  2178    385        4 (1653/1171-1912/1298) -21.0 W43-W69 S10-N00
      6 1786 1179 (S08W53 L=339)   341     61        0 (1773/1148-1800/1210)  36.9 W50-W54 S09-S06
      7  780 1107 (N16E12 L=275)   422    131        1 ( 756/1052- 813/1165) -81.9 E13-E09 N12-N19
      8 1830 1117 (S13W55 L=342)   494     80        1 (1814/1070-1845/1156)  45.0 W53-W57 S15-S10
      9  303  995 (N21E48 L=239)   289     73        0 ( 282/ 972- 333/1034) -59.0 E49-E45 N18-N23
     10 1692  775 (S28W36 L=322)  2414    543        0 (1552/ 551-1820/ 951)  32.2 W21-W50 S39-S21
     11  679  887 (N06E23 L=264)   637    117        2 ( 627/ 857- 737/ 926) -23.2 E26-E19 N04-N08
     12 1340  689 (S22W09 L=296)   529    154        0 (1287/ 634-1393/ 751)  26.6 W04-W13 S24-S19
     13  204  667 (N04E69 L=218)   576     72        1 ( 172/ 650- 244/ 688) -11.3 E73-E63 N02-N05
     14  694  585 (S11E31 L=256)  1287    331        0 ( 620/ 529- 774/ 679) -79.7 E34-E25 S15-S04
     15 1245  594 (S25W01 L=288)   265     71        2 (1227/ 552-1268/ 627)  78.7 W00-W02 S27-S22
     16  366  451 (S12E64 L=223)  4968    415        8 ( 258/ 306- 500/ 562) -81.5 E73-E56 S23-S03
     17  526  511 (S11E45 L=241)   300     42        0 ( 502/ 490- 551/ 532) -56.3 E46-E44 S12-S09
     18  989  308 (S36E26 L=261)   925    118        2 ( 942/ 269-1046/ 353)   0.0 E29-E20 S39-S33
  

    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