Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1391 1703 (N39W51 L=297)  1052    226        5 (1284/1671-1486/1723) -14.0 W40-W60 N35-N41
      2  590 1657 (N59E25 L=221)   942    196        3 ( 520/1597- 674/1714)  18.4 E33-E14 N54-N62
      3 1523 1667 (N32W59 L=306)   470     47        1 (1496/1646-1557/1683) -15.9 W56-W63 N31-N33
      4 1291 1477 (N28W29 L=275)   312     60        1 (1263/1454-1330/1499)   0.0 W25-W31 N26-N28
      5  356 1367 (N40E47 L=199)   800    132        0 ( 309/1342- 415/1391)  -8.1 E53-E39 N38-N42
      6 1765  835 (S22W47 L=294)   422     58        3 (1740/ 802-1790/ 865)  45.0 W44-W50 S23-S19
      7  485  565 (S11E46 L=200)   490     86        0 ( 440/ 540- 522/ 591)   7.1 E49-E41 S12-S09
      8  615  546 (S14E37 L=210)   400    115        0 ( 597/ 500- 642/ 587)  80.5 E38-E34 S17-S11
      9  583  453 (S20E43 L=203)   920    186        4 ( 513/ 418- 671/ 491) -15.3 E48-E35 S22-S16
     10  573  331 (S29E53 L=194)   469     87        1 ( 557/ 294- 590/ 370)  29.1 E57-E47 S31-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