Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  652 1527 (N33E22 L= 11)  4581    685       12 ( 401/1321- 890/1669) -31.1 E53-E05 N17-N46
      2 1355 1618 (N33W30 L= 63)   460    155        1 (1284/1595-1421/1645) -24.0 W25-W34 N31-N35
      3 1161 1488 (N26W13 L= 46)   772    310        1 (1060/1399-1291/1587) -45.0 W07-W20 N19-N33
      4 1679 1517 (N24W55 L= 87)   372     47        1 (1658/1496-1706/1533)   9.5 W52-W58 N22-N25
      5  345 1457 (N32E49 L=344)   415     65        0 ( 318/1432- 378/1480)  33.7 E52-E45 N30-N33
      6  435 1415 (N28E39 L=354)   451    130        1 ( 385/1368- 469/1457) -49.4 E44-E35 N24-N31
      7  394 1308 (N21E41 L=352)   752    205        1 ( 369/1254- 441/1399) -63.4 E43-E36 N17-N27
      8  240 1269 (N20E56 L=337)   992    150        2 ( 211/1199- 276/1331)  66.8 E59-E52 N16-N24
      9  183 1208 (N17E62 L=331)   838    138        2 ( 165/1135- 204/1260) -82.9 E65-E59 N12-N20
     10  479 1159 (N11E34 L=359)   575    255        1 ( 407/1075- 567/1229)  26.6 E39-E27 N06-N15
     11  311 1200 (N15E48 L=345)   423     69        1 ( 280/1175- 349/1223) -26.6 E50-E44 N13-N16
     12  329 1113 (N10E46 L=347)   598    122        2 ( 286/1092- 398/1151) -16.7 E49-E39 N08-N12
     13  146 1076 (N09E66 L=327)  2366    228        8 (  92/ 978- 230/1130)  33.7 E77-E55 N03-N13
     14  777  751 (S16E18 L= 15)   737    316        1 ( 686/ 647- 872/ 853)  33.7 E25-E10 S21-S10
     15  500  723 (S15E37 L=356)  1380    263        3 ( 474/ 592- 531/ 842)  72.6 E40-E33 S23-S07
     16 1469  730 (S23W27 L= 60)  1138    280        3 (1368/ 679-1583/ 823) -31.0 W19-W36 S26-S16
     17 1657  790 (S20W41 L= 74)   552    138        1 (1598/ 770-1732/ 807)   4.8 W36-W48 S21-S19
     18 1675  479 (S42W53 L= 86)   611     74        0 (1645/ 451-1709/ 515) -21.0 W47-W60 S43-S38
     19  356  449 (S31E64 L=329)   704     66        0 ( 339/ 408- 374/ 480)  29.1 E68-E59 S33-S28
  

    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