Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1027 1578 (N34W08 L=193)   235     82        2 (1011/1544-1042/1607)  78.7 W07-W09 N31-N36
      2 1259 1391 (N18W19 L=205)   195     71        0 (1240/1356-1290/1411)  45.0 W17-W21 N16-N19
      3 1397 1361 (N15W28 L=213)   205    113        0 (1346/1327-1444/1380)  14.0 W23-W31 N13-N15
      4 1777 1288 (N06W55 L=241)   556    130        0 (1756/1226-1808/1352)  82.9 W53-W57 N02-N10
      5  180 1005 (N09E60 L=125)   945    182        1 ( 150/ 939- 222/1080) -66.0 E63-E56 N04-N13
      6  282  899 (N01E51 L=134)   836    212        1 ( 242/ 810- 305/ 988)  48.4 E57-E48 S03-N06
      7 1321  917 (S10W16 L=202)   201     74        0 (1278/ 893-1347/ 935) -45.0 W13-W18 S11-S08
      8 1077  746 (S18E00 L=185)   643    200        1 ( 997/ 660-1168/ 815)  32.5 E06-W05 S21-S14
      9  280  646 (S14E58 L=127)  1300    267        3 ( 228/ 547- 315/ 731) -80.5 E61-E55 S19-S07
  

    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