Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  710 1702 (N48E32 L=342)   492    114        0 ( 688/1661- 728/1739) -54.5 E34-E28 N44-N51
      2  665 1615 (N41E32 L=343)   356     76        2 ( 642/1594- 690/1636)  33.7 E33-E29 N39-N42
      3 1035 1481 (N32E01 L= 14)   550    136        0 ( 977/1430-1090/1514) -40.6 E04-W03 N28-N34
      4 1309 1415 (N28W19 L= 33)   275     48        1 (1290/1386-1331/1441) -76.0 W17-W20 N26-N30
      5 1560 1323 (N22W37 L= 51)  1841    457        5 (1395/1237-1702/1414) -26.6 W25-W47 N16-N28
      6  600 1351 (N22E30 L=344)   420    102        1 ( 582/1309- 621/1411) -71.6 E31-E28 N19-N25
      7 1136 1276 (N19W07 L= 21)   849    196        5 (1082/1136-1192/1377) -71.6 W03-W10 N10-N25
      8 1298 1259 (N18W17 L= 32)   520    149        1 (1239/1194-1358/1328) -45.0 W13-W21 N14-N22
      9  728 1280 (N18E20 L=354)   222     73        0 ( 696/1261- 770/1302)  21.8 E22-E17 N16-N19
     10  942 1063 (N05E05 L=  9)  1295    352        1 ( 857/ 969-1089/1195) -40.2 E10-W03 N00-N13
     11  788  921 (S04E14 L=  0)   321    107        1 ( 738/ 899- 847/ 939)  15.9 E17-E10 S05-S02
     12 1113  245 (S53W12 L= 27)   355     75        0 (1071/ 231-1154/ 255) -11.3 W07-W17 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