Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  949 1707 (N53W06 L=187)   327     65        1 ( 923/1685- 976/1724)  21.8 W03-W09 N51-N54
      2 1377 1470 (N28W32 L=213)   240     75        0 (1360/1433-1400/1505) -90.0 W30-W33 N25-N30
      3  360 1415 (N37E51 L=130)   717     89        2 ( 326/1380- 396/1444) -16.7 E55-E45 N34-N39
      4  640 1322 (N29E23 L=158)   269     61        0 ( 610/1299- 667/1344)  21.8 E25-E20 N27-N30
      5  525 1299 (N28E32 L=149)   251     44        2 ( 500/1270- 549/1320)  56.3 E34-E30 N26-N29
      6  302 1151 (N20E52 L=130)   316     86        0 ( 274/1120- 334/1177) -33.7 E54-E48 N18-N22
      7  281 1064 (N15E53 L=128)   515    163        1 ( 256/ 992- 312/1148) -73.3 E56-E50 N10-N20
      8  169  926 (N07E68 L=113)   995     95        2 ( 142/ 884- 198/ 971)  31.0 E74-E63 N03-N09
      9 1446  359 (S47W26 L=208)  1288    207        2 (1342/ 329-1526/ 403)   0.0 W15-W34 S49-S43
     10  516  309 (S38E60 L=121)   951    108        1 ( 482/ 273- 550/ 355)  17.7 E75-E50 S41-S33
  

    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