Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1258 1762 (N56W09 L=167)   597     58        0 (1231/1733-1284/1797)  54.5 W06-W12 N52-N60
      2 1526 1708 (N56W43 L=201)   440     36        0 (1511/1683-1548/1728)  26.6 W38-W48 N53-N58
      3  355 1463 (N19E57 L=100)  2134    300        2 ( 284/1352- 447/1582)  76.0 E61-E53 N11-N27
      4  872 1528 (N30E18 L=140)   442     89        2 ( 844/1482- 896/1569) -56.3 E20-E16 N27-N33
      5  428 1483 (N21E51 L=107)   383     54        1 ( 404/1452- 452/1511)  78.7 E52-E49 N18-N23
      6  257 1461 (N17E71 L= 87)  1503    119        3 ( 221/1402- 294/1508) -74.1 E76-E66 N13-N20
      7 1566 1409 (N32W35 L=193)   833    198        3 (1476/1354-1645/1456) -12.1 W27-W42 N28-N34
      8  900  930 (S07E06 L=151)   348     70        1 ( 872/ 908- 934/ 956) -26.6 E08-E04 S08-S05
      9  200  895 (S18E62 L= 96)   844     67        1 ( 172/ 858- 229/ 932)  35.5 E66-E57 S20-S15
     10  555  716 (S25E28 L=130)   244     45        0 ( 531/ 698- 582/ 734) -26.6 E29-E25 S25-S23
     11 1241  517 (S28W22 L=180)   742    303        0 (1094/ 468-1344/ 576)  33.7 W13-W28 S33-S23
     12 1574  517 (S24W49 L=207)   472     64        2 (1548/ 482-1602/ 548) -78.7 W46-W50 S26-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