Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1227 1784 (N48E14 L=107)   386     60        1 (1198/1756-1249/1809) -26.6 E18-E12 N45-N49
      2  661 1700 (N23E46 L= 75)   371     43        0 ( 637/1681- 679/1725) -63.4 E47-E43 N21-N25
      3  553 1661 (N18E52 L= 69)   685    103        1 ( 517/1618- 579/1706) -74.1 E54-E49 N14-N21
      4 1546 1680 (N53W18 L=139)   444     51        0 (1518/1661-1579/1697)  32.0 W14-W22 N50-N56
      5  499 1648 (N16E55 L= 66)   558     63        1 ( 479/1618- 535/1674)  76.0 E58-E52 N14-N18
      6  627 1513 (N10E38 L= 83)   349     50        0 ( 603/1480- 648/1541)  90.0 E39-E36 N08-N12
      7 1349 1510 (N31W05 L=126)   228     50        0 (1328/1487-1365/1537)  76.0 W04-W06 N28-N32
      8 1830 1205 (N29W51 L=171)   425     96        0 (1790/1168-1866/1237)  -5.7 W45-W55 N28-N30
      9 1196  430 (S35W32 L=153)   590    124        1 (1163/ 369-1228/ 498) -32.5 W27-W38 S37-S30
     10 1359  392 (S31W46 L=166)   274     58        0 (1337/ 370-1386/ 410)  90.0 W44-W47 S32-S29
     11  835  369 (S52W12 L=133)   269     62        0 ( 809/ 350- 864/ 394)  -9.5 W08-W14 S52-S50
  

    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