Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1242 1808 (N51E03 L=286)   427     61        0 (1217/1792-1281/1824)  11.3 E05-W00 N49-N53
      2 1123 1628 (N32E06 L=283)   422    115        0 (1110/1560-1140/1683) -71.6 E08-E03 N27-N36
      3  524 1624 (N22E48 L=242)   355     73        0 ( 489/1604- 561/1638)  -8.1 E50-E43 N20-N22
      4 1751 1259 (N23W45 L=334)   593     97        2 (1719/1216-1784/1286)  56.3 W41-W48 N19-N25
      5 1029 1224 (N05E03 L=286)   278     73        1 (1000/1195-1069/1249)  53.1 E04-E01 N02-N06
      6  940  761 (S23E00 L=289)   259    102        0 ( 901/ 744- 990/ 786)  36.9 E02-W02 S24-S21
      7  989  554 (S35W08 L=297)   447    139        0 ( 929/ 538-1062/ 569)  32.0 W03-W12 S37-S32
      8  669  525 (S44E16 L=274)   576     93        1 ( 632/ 492- 706/ 560)  -8.1 E19-E12 S46-S41
      9 1197  236 (S52W43 L=333)   428     51        0 (1175/ 203-1220/ 261) -50.2 W40-W46 S55-S49
  

    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