Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1525 1528 (N29W47 L= 62)   792    230        0 (1456/1458-1603/1560)  -3.8 W37-W52 N25-N31
      2 1146 1442 (N30W15 L= 30)   634    230        2 (1122/1342-1172/1537)  77.0 W12-W17 N23-N36
      3 1625 1402 (N19W49 L= 65)   506    120        0 (1591/1353-1656/1460) -82.9 W48-W50 N15-N23
      4  420 1294 (N29E41 L=334)  1992    382        4 ( 384/1140- 491/1459) -62.4 E46-E34 N19-N40
      5 1522  984 (S04W31 L= 46)   424    141        0 (1494/ 921-1564/1037)  56.3 W28-W34 S07-S01
      6 1744  860 (S15W48 L= 64)   617     76        1 (1708/ 837-1772/ 889)   8.1 W44-W51 S17-S14
      7  715  576 (S18E26 L=349)   298     70        0 ( 683/ 553- 752/ 599)   9.5 E29-E23 S19-S17
      8  619  557 (S18E34 L=341)   380     71        1 ( 595/ 535- 654/ 587) -45.0 E36-E31 S20-S16
      9  563  510 (S21E40 L=335)   331     91        1 ( 538/ 478- 606/ 536)   0.0 E42-E36 S23-S19
     10  944  486 (S27E13 L=  2)   309    129        0 ( 885/ 474-1006/ 496)  -6.3 E17-E08 S28-S26
  

    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