Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  597 1626 (N17E45 L=277)   368     58        1 ( 578/1598- 622/1660) -63.4 E46-E42 N15-N19
      2  557 1569 (N13E45 L=277)   313     53        1 ( 527/1549- 584/1588)  33.7 E46-E42 N11-N14
      3 1325 1161 (N08W13 L=334)   504    123        0 (1272/1097-1359/1221) -32.0 W08-W16 N05-N10
      4  511 1122 (S14E33 L=288)   668    227        0 ( 410/1074- 593/1166)  32.5 E39-E27 S17-S10
      5 1222 1117 (N03W08 L=330)   288    110        0 (1176/1085-1283/1159)  -8.1 W05-W12 N02-N04
      6 1677 1003 (N10W39 L=  1)   317     65        0 (1662/ 965-1692/1039) -53.1 W37-W41 N08-N12
      7 1539  912 (N01W33 L=354)   806    269        0 (1464/ 802-1579/1016) -24.8 W25-W38 S03-N04
      8  831  745 (S28E04 L=318)   869    317        0 ( 684/ 707-1002/ 779)  16.7 E13-W07 S30-S24
      9 1599  327 (S26W66 L= 27)   601    164        1 (1568/ 310-1630/ 355)  23.2 W62-W70 S27-S24
  

    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