Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  570 1575 (N15E43 L=141)   289    124        0 ( 527/1554- 618/1599)  -7.1 E46-E38 N14-N16
      2  317 1465 (N05E59 L=125)  1412    206        6 ( 250/1415- 378/1535)  48.8 E64-E55 N01-N09
      3 1694 1287 (N26W36 L=220)   493     80        0 (1660/1269-1728/1311)  26.6 W33-W39 N23-N27
      4 1358 1108 (N05W17 L=201)   386    126        0 (1294/1080-1405/1143)  50.2 W13-W19 N02-N08
      5  865 1027 (S11E09 L=175)   218     78        1 ( 850/ 979- 882/1069) -59.0 E10-E07 S13-S08
      6  622  905 (S23E22 L=162)   669    254        0 ( 510/ 861- 736/ 970)  37.6 E29-E15 S26-S16
      7 1596  495 (S21W53 L=237)   368    107        2 (1563/ 457-1620/ 523)  90.0 W50-W54 S23-S18
      8 1406  485 (S27W40 L=224)   461    113        3 (1361/ 459-1458/ 511)  51.3 W37-W42 S29-S24
      9  863  203 (S66W24 L=208)   479     57        1 ( 836/ 184- 890/ 225)  45.0 W20-W27 S68-S62
  

    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