Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1496 1779 (N34W63 L= 42)   594    101        0 (1475/1743-1520/1812)  68.2 W61-W64 N31-N36
      2 1409 1565 (N26W42 L= 21)   620     84        1 (1369/1540-1452/1588) -26.6 W38-W44 N24-N28
      3 1774 1461 (N08W67 L= 47)  2609    378        0 (1673/1313-1821/1550) -90.0 W59-W74 N00-N15
      4 1440 1410 (N17W36 L= 16)   472     96        1 (1390/1382-1488/1431) -38.7 W33-W38 N14-N19
      5 1848 1131 (S11W59 L= 38)   767    105        3 (1828/1086-1875/1176)  51.3 W56-W61 S13-S08
      6  414 1068 (N24E38 L=301)   840    220        1 ( 352/ 986- 499/1144) -53.1 E42-E32 N17-N29
      7 1684 1077 (S08W42 L= 21)   498    148        0 (1666/1006-1715/1143)  49.4 W38-W45 S12-S05
      8  843 1039 (N12E10 L=330)   698    231        0 ( 731/ 971- 963/1095) -45.0 E15-E04 N05-N18
      9  261  822 (N12E57 L=282)   673    161        1 ( 213/ 801- 312/ 844) -24.0 E61-E52 N09-N13
     10 1000  722 (S09E09 L=330)   566    259        0 ( 878/ 671-1102/ 790) -47.7 E14-E04 S14-S03
     11  688  506 (S15E35 L=305)  2525    528        5 ( 632/ 333- 745/ 659)  53.1 E46-E28 S25-S05
     12  935  384 (S29E24 L=315)   357    102        0 ( 914/ 346- 957/ 423)  45.0 E27-E21 S31-S26
  

    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