Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  469 1732 (N30E68 L=345)   917    102        0 ( 447/1679- 483/1783) -24.0 E82-E63 N25-N33
      2 1592 1510 (N42W31 L= 85)   466     79        0 (1551/1477-1625/1535)  -8.1 W27-W34 N39-N42
      3 1126 1259 (N13W01 L= 54)   289     65        1 (1102/1230-1161/1285)   0.0 W00-W02 N11-N15
      4  186 1242 (S06E67 L=347)   506     44        2 ( 164/1220- 205/1268)  63.4 E70-E63 S07-S03
      5  777 1166 (N01E17 L= 36)  1089    257        0 ( 707/1065- 855/1268)  74.1 E19-E14 S06-N08
      6 1504 1211 (N19W26 L= 79)   329    103        0 (1450/1177-1544/1235)   0.0 W21-W28 N17-N19
      7  998  931 (S08W00 L= 54)   194     70        0 ( 970/ 902-1030/ 953)  76.0 W00-W01 S10-S06
      8 1582  799 (S03W40 L= 93)   689    111        1 (1539/ 766-1618/ 851)  59.0 W36-W41 S05-N00
      9 1862  774 (N02W68 L=121)  1929    172        3 (1804/ 719-1907/ 830) -10.0 W60-W77 N00-N05
  

    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