Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  613 1380 (N23E30 L= 75)   731    231        1 ( 505/1361- 719/1406)  -3.4 E39-E22 N21-N24
      2 1100 1206 (N14W04 L=109)   433    199        0 (1031/1178-1163/1252) -23.2 W00-W08 N12-N16
      3  314 1220 (N11E52 L= 53)   341     42        1 ( 299/1195- 331/1249)  90.0 E53-E50 N09-N13
      4 1316  744 (S14W20 L=125)   641    221        3 (1256/ 662-1384/ 819)  70.0 W15-W23 S19-S08
      5 1753  764 (S12W54 L=159)   390     70        1 (1723/ 737-1786/ 786)  26.6 W51-W57 S13-S10
      6  242  678 (S23E63 L= 42)  1038     75        3 ( 212/ 647- 277/ 714) -15.9 E67-E57 S25-S21
      7  683  551 (S29E23 L= 82)   241     59        2 ( 659/ 519- 704/ 572) -63.4 E24-E20 S31-S27
      8 1317  500 (S30W24 L=128)   386    107        0 (1293/ 450-1345/ 539)  71.6 W22-W24 S34-S27
      9  789  409 (S40E16 L= 89)   430     81        1 ( 768/ 374- 813/ 446) -90.0 E17-E13 S42-S36
     10 1394  266 (S50W43 L=148)   783     83        4 (1354/ 234-1429/ 297) -17.4 W36-W52 S53-S47
  

    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