Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  441 1465 (N46E32 L=116)   436     43        0 ( 419/1436- 470/1487)  33.7 E34-E28 N44-N48
      2 1044 1339 (N22W10 L=158)   218     84        0 (1001/1315-1082/1352)   0.0 W06-W11 N20-N22
      3  368 1120 (N26E41 L=107)   840    250        1 ( 315/1048- 449/1167)   4.4 E46-E33 N22-N28
      4 1705  913 (S20W40 L=188)   648    169        3 (1668/ 851-1772/ 966)  16.7 W36-W46 S22-S18
      5 1632  687 (S32W31 L=179)   517    185        0 (1553/ 641-1701/ 738)   8.1 W23-W37 S32-S30
      6 1260  600 (S26W03 L=151)   208     82        0 (1232/ 560-1283/ 632)  36.9 W00-W04 S27-S24
      7  859  572 (S17E22 L=126)   251     89        0 ( 816/ 543- 894/ 606)  15.9 E25-E18 S17-S15
      8  450  524 (S10E51 L= 97)   835    147        1 ( 425/ 451- 478/ 586)  37.9 E55-E46 S14-S06
      9 1426  499 (S38W11 L=159)   589    173        0 (1358/ 461-1499/ 553)   8.7 W05-W18 S38-S36
     10  905  311 (S35E31 L=117)   318     57        0 ( 886/ 283- 922/ 344)  23.2 E34-E27 S36-S32
     11  560  287 (S28E60 L= 88)  1704    141        2 ( 504/ 244- 618/ 341) -69.4 E64-E54 S31-S23
  

    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