Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1385 1526 (N34W33 L=355)   771    183        1 (1311/1485-1481/1566) -32.5 W27-W39 N29-N37
      2  794 1540 (N40E14 L=309)   298     60        0 ( 768/1517- 820/1562)  26.6 E15-E11 N38-N41
      3 1485 1339 (N20W34 L=357)   681    184        2 (1440/1254-1539/1405)  36.0 W29-W40 N15-N23
      4 1446 1191 (N12W29 L=351)   243     73        0 (1423/1160-1474/1221) -63.4 W27-W30 N09-N13
      5 1854 1198 (N07W66 L= 28)   643     60        2 (1826/1171-1881/1221)   5.7 W61-W71 N05-N07
      6 1667  852 (S12W42 L=  5)   703    208        1 (1635/ 754-1691/ 953)  76.0 W39-W44 S17-S05
      7  463  554 (S22E45 L=277)   504     61        0 ( 432/ 526- 487/ 574) -33.7 E47-E42 S23-S20
      8  944  500 (S28E10 L=313)   376    125        0 ( 914/ 442- 961/ 546) -82.9 E11-E08 S32-S24
      9  749  437 (S31E26 L=297)   741    185        0 ( 722/ 351- 795/ 510)  50.2 E30-E20 S38-S26
     10  680  294 (S43E39 L=283)   570     85        0 ( 646/ 262- 712/ 324) -71.6 E41-E37 S46-S40
  

    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