Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  664 1400 (N33E14 L=179)  1039    206        2 ( 590/1332- 732/1444) -37.9 E19-E10 N28-N36
      2 1545 1326 (N08W39 L=232)   550    120        1 (1514/1289-1589/1356) -21.8 W36-W41 N05-N09
      3  154 1019 (N20E62 L=131)  1758    297        4 ( 114/ 893- 180/1121) -82.4 E66-E58 N12-N27
      4  851  784 (S07E15 L=178)  3474    731        9 ( 671/ 560- 995/ 925) -10.0 E26-E07 S20-N03
      5 1252  868 (S11W09 L=203)   424    169        0 (1189/ 852-1331/ 885) -20.6 W05-W13 S12-S09
      6 1277  762 (S17W09 L=202)   331    143        1 (1220/ 732-1354/ 789)  -6.3 W04-W13 S18-S16
      7  674  640 (S11E30 L=163)   364    119        1 ( 642/ 592- 708/ 690)  29.7 E33-E26 S13-S08
      8 1059  558 (S25E09 L=184)   506    138        0 ( 998/ 524-1139/ 602)  11.3 E13-E03 S25-S23
      9  814  530 (S21E25 L=168)   454    107        1 ( 766/ 500- 865/ 561)   7.1 E28-E20 S21-S19
     10  514  306 (S29E58 L=135)   564     93        0 ( 489/ 271- 535/ 340)  71.6 E59-E55 S31-S25
  

    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