Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  529 1778 (N59E47 L=203)   933    126        0 ( 459/1764- 586/1793)  -2.1 E64-E37 N58-N60
      2 1669 1631 (N30W63 L=313)   829    101        3 (1631/1606-1709/1671) -39.8 W58-W70 N28-N33
      3 1345 1564 (N29W29 L=279)   851    168        1 (1278/1489-1403/1621) -73.3 W25-W31 N23-N33
      4  959 1458 (N27W01 L=250)   336     83        2 ( 921/1434-1008/1488) -36.9 E01-W03 N24-N28
      5 1674 1456 (N18W52 L=301)   450     80        0 (1635/1428-1709/1485) -38.7 W48-W54 N16-N20
      6 1629  909 (S14W38 L=288)   283     84        0 (1599/ 883-1646/ 948)  45.0 W35-W39 S15-S11
      7 1778  621 (S34W57 L=307)   438     44        0 (1763/ 595-1797/ 649) -45.0 W54-W61 S35-S31
      8 1538  584 (S34W32 L=282)   427     61        1 (1510/ 551-1554/ 623) -78.7 W29-W33 S36-S31
      9  961  472 (S34E12 L=238)   492    119        1 ( 934/ 422- 992/ 528) -76.0 E12-E09 S38-S30
     10  967  248 (S53E19 L=230)   450     81        1 ( 942/ 219-1018/ 270)  16.7 E23-E13 S54-S50
     11  806  170 (S57E45 L=204)   902     77        2 ( 774/ 141- 833/ 207)  13.2 E55-E38 S60-S53
  

    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