Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  546  981 (N16E30 L=  5)  2215    539        5 ( 440/ 842- 715/1179) -64.4 E34-E21 N04-N29
      2  893 1104 (N15E05 L= 29)   647    259        4 ( 784/1062- 984/1139) -38.7 E10-W00 N10-N19
      3  274  743 (N07E58 L=336)   535     91        0 ( 239/ 713- 301/ 787)  23.2 E62-E54 N05-N09
      4 1644  750 (S27W32 L= 66)   514     96        0 (1606/ 719-1682/ 781)   8.1 W28-W35 S28-S25
      5 1404  510 (S35W09 L= 44)  1374    326        0 (1271/ 429-1477/ 616)  22.4 E01-W16 S38-S29
      6  681  535 (S13E34 L=  0)   372    103        0 ( 647/ 489- 711/ 587)  81.9 E35-E33 S15-S08
      7 1139  271 (S44E18 L= 16)   577     73        1 (1101/ 245-1175/ 302)  15.9 E21-E14 S46-S42
      8 1089  234 (S46E25 L= 10)   444     36        1 (1068/ 215-1114/ 259) -90.0 E26-E22 S48-S42
  

    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