Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  376 1659 (N37E56 L=230)   718    108        1 ( 355/1617- 397/1693) -80.5 E58-E54 N33-N39
      2  989 1427 (N20E03 L=283)   454     69        3 ( 951/1400-1030/1447) -11.3 E05-W00 N18-N21
      3 1434 1363 (N18W25 L=311)   377     55        0 (1412/1320-1453/1398) -78.7 W23-W26 N15-N20
      4 1254 1285 (N12W13 L=299)   952    200        0 (1154/1218-1351/1363) -33.7 W07-W19 N08-N16
      5  671 1237 (N08E22 L=264)   305     95        0 ( 650/1189- 694/1287) -63.4 E23-E20 N05-N11
      6  439 1216 (N07E38 L=248)   517    138        1 ( 393/1171- 495/1271) -29.1 E42-E33 N03-N09
      7 1732 1230 (N12W47 L=333)   425     85        1 (1701/1199-1765/1268) -31.0 W44-W50 N09-N13
      8  999 1214 (N07E02 L=284)   342    143        0 ( 934/1191-1063/1249)   9.5 E05-W01 N05-N09
      9 1871 1015 (S00W61 L=348)   746    111        1 (1834/ 982-1900/1044) -14.0 W57-W65 S01-N01
     10  179  906 (S11E62 L=224)  1242    204        3 ( 137/ 863- 241/ 962)  14.0 E67-E54 S13-S07
     11 1863  881 (S08W62 L=348)  1642    270        3 (1814/ 817-1901/ 953)  56.3 W56-W67 S12-S03
     12  367  698 (S25E46 L=240)  2638    580        4 ( 322/ 458- 460/ 875)  56.3 E56-E36 S40-S13
     13 1512  819 (S14W32 L=318)   338    113        1 (1481/ 763-1561/ 865)  63.4 W30-W35 S17-S11
     14  654  846 (S16E23 L=263)   199     57        0 ( 618/ 835- 685/ 861)   0.0 E25-E20 S16-S14
     15 1192  632 (S27W13 L=299)   710    140        2 (1115/ 593-1252/ 687)  29.1 W07-W16 S30-S23
     16 1530  609 (S27W37 L=323)   497    111        2 (1477/ 568-1565/ 632)   0.0 W32-W39 S29-S25
     17 1699  482 (S34W56 L=342)   407    106        0 (1675/ 470-1733/ 500)   0.0 W53-W59 S34-S32
  

    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