Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1022 1615 (N31E17 L=232)  2024    342        5 ( 911/1517-1113/1685) -14.7 E27-E08 N25-N34
      2  806 1569 (N22E28 L=221)  1295    242        3 ( 741/1486- 894/1650)  77.9 E30-E26 N15-N29
      3 1518 1311 (N26W23 L=273)  2707    548        4 (1381/1074-1586/1485) -52.7 W10-W31 N12-N35
      4  172 1140 (S16E65 L=184)  3577    318        6 (  84/1016- 244/1272)  38.7 E84-E56 S24-S07
      5  411  957 (S22E37 L=212)   849    267        0 ( 355/ 872- 516/1008)  47.7 E41-E30 S27-S16
      6 1788  958 (N13W52 L=301)   451     87        0 (1772/ 911-1802/1007)  90.0 W50-W52 N09-N15
      7  871  843 (S18E04 L=245)   368    178        0 ( 835/ 780- 921/ 883)  90.0 E06-E01 S21-S14
      8  436  850 (S28E34 L=215)   495     96        1 ( 390/ 835- 492/ 868)  14.0 E37-E29 S29-S26
      9 1335  775 (S10W24 L=274)   392    169        0 (1272/ 719-1393/ 802)  -6.3 W20-W29 S11-S08
     10  756  383 (S51W01 L=251)  1609    238        3 ( 659/ 309- 841/ 472) -16.7 E09-W11 S54-S46
     11  956  241 (S55W28 L=277)   335     36        0 ( 934/ 222- 976/ 258)  45.0 W25-W30 S56-S53
  

    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