Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  733 1848 (N40E47 L= 22)  1324    235        2 ( 687/1809- 769/1883) -18.4 E53-E41 N36-N43
      2  779 1675 (N27E32 L= 37)   551    183        2 ( 752/1596- 804/1731) -48.4 E36-E28 N21-N31
      3  527 1444 (N08E40 L= 29)   595    220        2 ( 437/1400- 609/1478)  -8.1 E47-E33 N06-N09
      4 1260 1370 (N18W06 L= 75)   432    220        1 (1210/1270-1298/1458) -51.3 W01-W10 N12-N22
      5  513 1200 (S07E35 L= 34)  1144    318        1 ( 420/1051- 601/1358)  75.3 E38-E32 S16-N03
      6 1580 1186 (N15W30 L= 99)   478    186        0 (1553/1089-1628/1259) -56.3 W27-W35 N10-N19
      7 1151 1037 (S04W07 L= 76)   331    135        1 (1101/ 979-1206/1105)  69.4 W04-W08 S07-N01
      8  287  965 (S24E49 L= 20)   475    115        1 ( 245/ 945- 322/ 988)  15.9 E52-E45 S25-S21
  

    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