Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  388 1416 (N43E42 L=267)   560    126        2 ( 349/1385- 429/1443) -29.1 E46-E37 N40-N45
      2  287 1360 (N40E51 L=258)   832    148        1 ( 266/1313- 328/1415) -60.3 E53-E47 N37-N44
      3 1293 1346 (N20W24 L=333)   950    250        2 (1205/1298-1399/1376) -32.5 W18-W29 N15-N23
      4 1657 1331 (N10W49 L=358)   353     66        0 (1626/1308-1688/1353) -36.9 W46-W51 N07-N11
      5  748 1309 (N29E12 L=297)   521    114        3 ( 713/1277- 801/1342)   8.1 E14-E07 N27-N31
      6  280 1169 (N28E50 L=258)   996    178        2 ( 255/1091- 311/1241) -66.0 E52-E47 N23-N32
      7 1394 1040 (N00W22 L=331)   249     95        1 (1347/ 999-1428/1076)  15.9 W18-W25 S01-N01
      8 1387  893 (S08W19 L=328)   722    244        1 (1309/ 838-1496/ 971) -38.7 W15-W25 S12-S02
      9  653  894 (N06E25 L=284)   294    106        2 ( 592/ 882- 699/ 910)  -8.1 E28-E21 N04-N07
     10  920  555 (S19E16 L=293)  1446    405        5 ( 800/ 399- 989/ 686)  38.3 E28-E09 S27-S11
     11  374  550 (S11E56 L=253)  1241    326        3 ( 307/ 514- 441/ 607)   7.1 E64-E48 S13-S07
     12  455  391 (S23E57 L=252)   634    117        1 ( 422/ 366- 496/ 412)   7.1 E60-E51 S24-S21
     13  927  300 (S38E26 L=283)   845    155        2 ( 890/ 243- 972/ 338)  28.6 E31-E20 S42-S35
  

    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