Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  364 1520 (N34E47 L=327)   454     55        0 ( 344/1487- 389/1548)  63.4 E49-E45 N31-N36
      2 1877 1291 (N10W66 L= 80)  1408    283        5 (1847/1216-1917/1375) -68.2 W63-W71 N05-N15
      3  304 1255 (N17E48 L=326)   399     53        1 ( 283/1223- 319/1289)  90.0 E50-E46 N14-N18
      4 1050 1235 (N10W03 L= 17)   230     50        1 (1030/1210-1065/1270)  56.3 W01-W03 N08-N12
      5 1445 1240 (N09W27 L= 41)   236     75        1 (1405/1228-1482/1256) -21.8 W24-W29 N07-N09
      6  502 1185 (N11E32 L=342)   348     95        1 ( 460/1165- 551/1205)   0.0 E35-E28 N09-N12
      7  626 1087 (N04E24 L=350)   750    200        0 ( 567/ 980- 653/1161) -68.2 E27-E22 S02-N08
      8  276 1027 (N03E51 L=323)  1049    202        3 ( 244/ 962- 328/1084)  63.4 E53-E46 N00-N06
      9  365 1026 (N02E43 L=331)   435     87        2 ( 328/ 994- 408/1055) -18.4 E45-E39 N00-N04
     10  474  699 (S18E39 L=335)   241     56        1 ( 445/ 683- 502/ 717) -14.0 E41-E36 S18-S16
     11 1639  681 (S26W41 L= 55)   337     67        1 (1611/ 640-1671/ 711)  45.0 W39-W43 S28-S24
     12 1177  447 (S40W07 L= 21)   250     47        0 (1147/ 432-1208/ 462)  26.6 W05-W09 S40-S38
  

    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