Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  978 1588 (N40E05 L=342)   496    113        3 ( 951/1543-1001/1631) -74.1 E06-E02 N36-N43
      2 1543 1430 (N29W38 L= 24)   601    108        2 (1500/1406-1591/1451)   0.0 W33-W41 N27-N30
      3  350 1409 (N25E52 L=294)   374    126        0 ( 324/1371- 379/1447) -35.5 E55-E48 N22-N27
      4 1795 1360 (N24W60 L= 46)  2363    444        5 (1664/1274-1906/1425) -21.0 W47-W73 N17-N28
      5 1006 1324 (N22E02 L=345)   411     81        0 ( 973/1284-1030/1377) -56.3 E04-W00 N19-N25
      6 1745 1172 (N12W51 L= 38)   872    118        1 (1692/1103-1788/1247) -51.3 W47-W55 N07-N17
      7 1355 1064 (N06W21 L=  7)   219     63        0 (1340/1027-1371/1098) -78.7 W19-W21 N03-N08
      8 1491  983 (N01W30 L= 16)  1104    205        2 (1402/ 945-1562/1044)  18.4 W23-W35 S01-N04
      9 1256  902 (S04W15 L=  1)   360    179        0 (1167/ 879-1332/ 924)  -5.2 W08-W19 S05-S02
  

    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