Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  513 1659 (N52E47 L= 74)   474     52        1 ( 486/1633- 541/1682) -18.4 E53-E41 N50-N54
      2 1751 1441 (N22W64 L=185)   440     73        2 (1734/1406-1770/1472)  24.0 W59-W68 N19-N23
      3 1359 1108 (N07W22 L=142)   514    145        1 (1330/1045-1385/1171) -90.0 W20-W23 N03-N11
      4 1239  858 (S07W12 L=132)   206    107        0 (1188/ 827-1268/ 891)  31.0 W08-W13 S08-S04
      5 1352  773 (S13W18 L=139)   361    110        0 (1295/ 750-1403/ 798) -23.2 W14-W21 S15-S11
      6 1820  721 (S24W61 L=181)   932    130        2 (1771/ 668-1866/ 763)  20.0 W54-W67 S26-S21
      7  551  523 (S23E39 L= 82)   308     42        1 ( 524/ 504- 580/ 541) -33.7 E41-E36 S24-S21
      8  474  470 (S27E48 L= 72)   538     66        1 ( 455/ 430- 501/ 507)  35.5 E52-E44 S29-S24
      9 1617  452 (S41W45 L=166)   592     74        1 (1593/ 424-1634/ 496) -71.6 W41-W48 S43-S36
     10 1513  363 (S47W38 L=158)   413     56        0 (1497/ 335-1532/ 386) -59.0 W35-W41 S49-S44
  

    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