Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  809 1576 (N25E28 L=329)   648    128        0 ( 752/1534- 858/1600)  -6.3 E32-E23 N23-N27
      2 1685 1495 (N45W42 L= 39)   376     32        0 (1660/1482-1711/1508)   0.0 W38-W45 N43-N45
      3 1426 1425 (N31W17 L= 14)   434    170        0 (1396/1343-1460/1492) -52.1 W13-W20 N26-N35
      4  744 1305 (N07E23 L=334)   622    191        0 ( 692/1263- 800/1373) -26.6 E28-E19 N05-N10
      5 1077 1110 (N03W01 L=358)   225     82        0 (1042/1094-1119/1131)   0.0 E01-W03 N02-N04
      6 1155  681 (S19W16 L= 13)   469    167        2 (1081/ 643-1223/ 711)  45.0 W12-W19 S23-S16
      7 1333  477 (S27W36 L= 33)   375     76        0 (1291/ 456-1372/ 492)  45.0 W33-W38 S29-S25
      8 1385  346 (S34W49 L= 46)  1016    123        2 (1332/ 313-1444/ 389)   0.0 W42-W56 S35-S32
      9  757  282 (S60W03 L=360)   343     49        1 ( 733/ 259- 773/ 306) -45.0 W00-W05 S62-S57
  

    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