Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1522 1680 (N35W47 L=180)   343     65        0 (1498/1659-1542/1701)  21.8 W44-W49 N33-N36
      2  518 1503 (N35E31 L=101)  1948    543       10 ( 348/1428- 699/1577) -20.6 E48-E16 N27-N42
      3  268 1193 (N17E50 L= 82)   939    186        2 ( 236/1127- 316/1251)  54.5 E53-E45 N13-N20
      4  879 1174 (N09E07 L=125)   196     81        1 ( 841/1140- 916/1206) -45.0 E09-E05 N07-N11
      5 1081 1069 (N01W04 L=136)   235     64        0 (1060/1035-1103/1120)  45.0 W02-W05 N00-N04
      6 1795  872 (S18W52 L=185)   663    207        0 (1733/ 820-1858/ 938)  22.6 W46-W58 S20-S13
      7  939  789 (S14E08 L=125)   181    110        0 ( 898/ 741- 992/ 817)  33.7 E10-E04 S16-S12
      8  854  751 (S15E13 L=119)   193     75        1 ( 834/ 718- 887/ 785) -63.4 E14-E11 S17-S13
      9 1294  459 (S40W13 L=146)   738    175        2 (1207/ 412-1380/ 528)  17.1 W06-W19 S42-S35
     10  897  441 (S36E17 L=116)   508    153        0 ( 856/ 392- 926/ 511)  41.2 E21-E13 S38-S31
     11  589  371 (S37E43 L= 90)   989    327        0 ( 561/ 252- 628/ 477)  62.1 E48-E37 S46-S29
  

    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