Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  596 1545 (N30E37 L= 74)   478    101        3 ( 563/1517- 645/1569)  18.4 E39-E33 N28-N31
      2  742 1443 (N25E23 L= 88)  1593    389        5 ( 633/1325- 818/1517) -30.5 E33-E16 N17-N29
      3 1405 1474 (N34W23 L=135)   381     71        0 (1374/1436-1435/1515) -59.0 W21-W25 N31-N36
      4 1853 1178 (N18W62 L=173)   454     66        1 (1831/1160-1875/1199) -21.8 W59-W64 N16-N18
      5  715 1077 (N02E20 L= 92)  1084    285        2 ( 594/ 992- 825/1119) -10.6 E27-E11 S02-N04
      6  550 1034 (S02E30 L= 82)   590    124        1 ( 500/ 967- 591/1084)  45.0 E33-E27 S06-N00
      7  246  736 (S24E56 L= 55)   593    148        0 ( 223/ 683- 261/ 787)  58.0 E59-E54 S28-S20
      8  969  686 (S20E00 L=111)   197     63        0 ( 945/ 654- 999/ 715)  68.2 E01-W01 S22-S17
  

    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