Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1403 1741 (N53W43 L=207)   518     70        0 (1369/1716-1431/1768)  14.0 W36-W49 N50-N55
      2  988 1733 (N53E03 L=161)   869    185        1 ( 905/1706-1070/1759)  18.4 E12-W05 N50-N55
      3  376 1344 (N23E49 L=115)   314     85        0 ( 361/1307- 391/1386) -59.0 E51-E46 N20-N25
      4  489 1263 (N18E37 L=127)   428    146        0 ( 444/1229- 553/1290)  20.6 E40-E32 N15-N19
      5  225 1066 (N05E59 L=105)   505     69        0 ( 200/1036- 256/1097)  33.7 E62-E56 N02-N06
      6  185  970 (S01E65 L=100)   713     59        0 ( 153/ 948- 219/ 992)  -6.3 E69-E60 S02-N00
      7  534  571 (S26E36 L=128)  1497    427        1 ( 411/ 420- 662/ 710)  33.7 E57-E23 S38-S16
      8 1411  683 (S18W26 L=190)   258    104        0 (1355/ 667-1460/ 699)   7.1 W21-W29 S19-S17
      9  813  486 (S32E16 L=148)  3205    607        8 ( 606/ 349-1019/ 681)  34.8 E37-W00 S43-S17
     10 1620  515 (S31W48 L=212)   542     59        0 (1589/ 488-1648/ 546)  53.1 W45-W50 S32-S28
     11 1533  322 (S47W53 L=217)  1389     94        1 (1504/ 287-1565/ 366) -22.2 W45-W67 S51-S42
  

    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