Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  744 1569 (N36E14 L= 59)   548    107        1 ( 690/1545- 791/1588)   0.0 E18-E10 N34-N37
      2  268 1363 (N28E55 L= 18)  3074    380        6 ( 229/1204- 331/1550) -39.6 E78-E45 N17-N41
      3 1373 1459 (N22W28 L=101)   863    232        0 (1260/1429-1489/1481) -18.4 W20-W35 N19-N24
      4  569 1378 (N25E27 L= 46)   240     46        0 ( 549/1349- 586/1406) -63.4 E28-E25 N23-N27
      5 1856 1282 (N07W65 L=138)  2869    209        7 (1783/1221-1935/1359) -33.7 W55-W77 N02-N12
      6 1749  911 (S14W49 L=122)   309     27        0 (1731/ 889-1764/ 936)  56.3 W47-W50 S15-S12
      7  289  659 (S15E58 L= 15)   979    189        2 ( 239/ 581- 327/ 743)  29.4 E67-E51 S19-S10
      8  818  663 (S21E17 L= 56)   248     59        0 ( 788/ 631- 839/ 698) -63.4 E18-E15 S22-S18
      9  943  618 (S25E10 L= 64)   322     61        1 ( 911/ 602- 985/ 640) -21.8 E11-E06 S26-S23
     10 1333  447 (S41W17 L= 90)   636     79        2 (1286/ 424-1382/ 477) -14.0 W13-W21 S42-S38
     11 1223  415 (S42W08 L= 81)   460     61        1 (1203/ 387-1248/ 452)  78.7 W06-W10 S44-S39
     12 1136  380 (S44W01 L= 74)   420     67        1 (1103/ 364-1175/ 408) -11.3 E01-W04 S45-S41
  

    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