Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  583 1654 (N27E49 L=253)   327     84        1 ( 543/1634- 615/1674)  18.4 E51-E45 N25-N28
      2  825 1517 (N23E24 L=277)   373    148        0 ( 786/1455- 850/1564) -40.6 E28-E20 N19-N26
      3  534 1528 (N18E45 L=256)   310     57        2 ( 521/1492- 553/1554) -63.4 E47-E43 N15-N19
      4 1185 1454 (N27W00 L=302)   921    284        0 (1064/1364-1295/1506) -10.0 E08-W09 N23-N29
      5 1367 1427 (N30W13 L=315)   312     98        0 (1331/1388-1400/1470) -33.7 W10-W16 N27-N32
      6 1866 1184 (N27W62 L=  4)  2299    226        3 (1822/1084-1897/1294)  81.3 W56-W68 N20-N34
      7  425 1069 (S11E39 L=262)   623    104        2 ( 379/1038- 467/1095)  33.7 E42-E36 S13-S09
      8 1745 1001 (N12W49 L=350)   650    104        1 (1724/ 952-1773/1039)  56.3 W46-W51 N09-N15
      9 1319  719 (S13W25 L=326)   224    111        0 (1275/ 708-1374/ 740)   8.1 W21-W28 S14-S12
     10 1430  568 (S20W38 L=339)  1971    356        4 (1345/ 442-1569/ 670)  67.6 W33-W43 S28-S10
     11  454  558 (S44E35 L=267)   346     53        0 ( 440/ 528- 470/ 585)  78.7 E36-E32 S46-S41
     12  909  540 (S34W04 L=306)   671    205        1 ( 817/ 511- 997/ 567)  28.6 E02-W09 S37-S29
  

    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