Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1274 1545 (N35W02 L= 77)  1165    240        3 (1217/1440-1313/1642) -50.2 E04-W07 N28-N40
      2  400 1542 (N13E58 L= 17)  1459    319        2 ( 339/1449- 503/1641)  74.1 E62-E54 N06-N20
      3  655 1464 (N14E34 L= 41)   343    130        0 ( 636/1408- 694/1505) -81.9 E35-E32 N09-N16
      4 1581 1176 (N19W31 L=106)  1143    228        3 (1535/1071-1631/1274)  85.9 W28-W33 N12-N26
      5 1808 1037 (N18W53 L=128)  1581    196        3 (1767/ 948-1840/1130) -52.1 W48-W58 N12-N22
      6 1237 1040 (N03W12 L= 86)   389    127        1 (1169/1022-1293/1063)  29.7 W07-W14 N00-N05
      7 1895  873 (N11W70 L=144)   938     94        1 (1866/ 842-1920/ 928)  80.5 W65-W74 N08-N14
      8  917  746 (S22W01 L= 75)   532    229        2 ( 821/ 708-1017/ 779)  35.0 E05-W05 S26-S17
      9  525  713 (S34E26 L= 49)  1024    154        0 ( 433/ 673- 614/ 753)  14.9 E34-E19 S36-S29
     10 1540  558 (S17W47 L=122)   674    194        3 (1473/ 528-1606/ 602)  45.0 W41-W50 S19-S12
     11 1649  558 (S13W56 L=131)   380     59        0 (1618/ 539-1677/ 583)  45.0 W53-W58 S14-S11
     12 1464  427 (S26W49 L=124)   867    214        3 (1386/ 364-1539/ 485)  -3.2 W40-W58 S28-S23
  

    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