Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  859 1823 (N63E23 L= 14)  1144    174        1 ( 792/1803- 933/1848)  -2.9 E33-E12 N60-N65
      2 1155 1725 (N52W13 L= 50)   332     24        1 (1132/1704-1178/1744) -11.3 W10-W15 N50-N54
      3  560 1469 (N31E36 L=  1)   461    106        3 ( 521/1430- 598/1517) -29.1 E40-E31 N28-N34
      4 1200 1392 (N27W12 L= 49)   701    130        1 (1136/1352-1247/1425) -20.6 W07-W15 N24-N29
      5  865 1306 (N21E11 L= 26)   356    149        1 ( 817/1278- 912/1342) -23.2 E14-E07 N19-N23
      6 1750  853 (S08W52 L= 89)   519     36        1 (1726/ 828-1774/ 881) -21.8 W49-W54 S09-S06
      7 1407  758 (S13W25 L= 62)   635    189        0 (1379/ 699-1456/ 834) -71.6 W22-W29 S17-S08
      8 1491  750 (S14W31 L= 68)   234     62        0 (1470/ 716-1511/ 778)  63.4 W29-W32 S16-S12
      9  837  684 (S18E12 L= 25)   857    175        3 ( 766/ 645- 922/ 721)  20.0 E17-E06 S20-S15
     10  389  523 (S31E52 L=344)  2635    332        4 ( 314/ 435- 487/ 636) -57.3 E58-E45 S37-S22
     11 1413  569 (S26W28 L= 65)  1116    280        0 (1301/ 534-1541/ 602)   3.0 W19-W38 S28-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