Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  719 1701 (N48E31 L=342)   499    118        1 ( 696/1662- 737/1743) -54.5 E34-E27 N44-N51
      2  676 1619 (N41E31 L=342)   424     62        2 ( 651/1596- 701/1646)  45.0 E32-E28 N39-N43
      3 1054 1481 (N32W01 L= 14)   502    116        1 ( 998/1442-1103/1511) -29.7 E03-W04 N29-N34
      4 1320 1419 (N29W20 L= 33)   282     71        1 (1297/1385-1343/1444)  63.4 W18-W21 N26-N30
      5 1569 1324 (N22W38 L= 51)  1900    393        4 (1409/1238-1714/1414) -26.6 W26-W48 N16-N28
      6  615 1345 (N22E29 L=344)   422     60        1 ( 596/1308- 637/1382)  90.0 E30-E27 N19-N24
      7 1156 1269 (N19W08 L= 21)  1022    239        3 (1096/1129-1208/1375) -69.4 W03-W11 N09-N25
      8 1309 1267 (N19W18 L= 31)   723    189        0 (1244/1189-1386/1351) -51.3 W13-W23 N13-N23
      9  966 1059 (N05E04 L= 10)  1302    361        2 ( 872/ 968-1119/1198) -35.5 E09-W05 N00-N13
     10  813  921 (S04E13 L=  0)   306    109        1 ( 762/ 906- 868/ 936)   0.0 E16-E09 S04-S02
     11 1135  244 (S53W14 L= 28)   383     79        0 (1087/ 230-1171/ 255)  -6.3 W09-W18 S54-S51
  

    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