Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1098 1523 (N30E07 L=295)   366     82        1 (1063/1482-1119/1560) -80.5 E09-E05 N26-N32
      2 1595 1316 (N28W32 L=334)   936    137        2 (1537/1268-1653/1349)   0.0 W26-W37 N25-N30
      3  206  945 (S23E60 L=243)   875    195        2 ( 148/ 888- 274/1015)  33.7 E69-E51 S27-S17
      4 1618  967 (N07W38 L=341)   299     47        1 (1595/ 945-1639/ 992)  63.4 W36-W39 N05-N09
      5  476  911 (S20E33 L=269)   389     93        1 ( 435/ 880- 523/ 944)  45.0 E35-E29 S22-S17
      6 1794  824 (N03W58 L=  0)   771    104        3 (1761/ 767-1829/ 858)  45.0 W55-W61 N00-N05
      7  896  760 (S20E02 L=300)   413    144        2 ( 829/ 740- 967/ 778)  14.0 E06-W02 S22-S18
      8  265  728 (S36E54 L=248)   861     89        2 ( 222/ 702- 310/ 764) -20.6 E60-E48 S37-S34
      9 1428  383 (S31W48 L=351)   343    107        0 (1382/ 352-1471/ 415)  63.4 W46-W49 S34-S27
  

    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