Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1073 1613 (N43W12 L=131)   263     24        1 (1048/1598-1099/1630)   0.0 W09-W14 N42-N44
      2 1504 1588 (N35W47 L=166)   822    112        2 (1475/1536-1540/1625)  20.0 W42-W53 N31-N37
      3 1230 1383 (N26W18 L=137)   296    127        0 (1177/1344-1283/1427) -50.2 W15-W20 N22-N29
      4 1642 1308 (N15W47 L=166)  1385    261        7 (1581/1218-1697/1413) -85.6 W42-W51 N09-N22
      5  215 1277 (N26E66 L= 53)   914    131        0 ( 173/1254- 267/1310)   0.0 E75-E58 N24-N28
      6 1358 1239 (N15W24 L=143)   228     78        0 (1333/1213-1381/1278) -76.0 W22-W25 N13-N17
      7  878 1090 (N11E08 L=111)   597    228        2 ( 798/1052-1002/1124) -22.6 E13-E01 N07-N13
      8 1757  860 (S14W50 L=169)   447     52        1 (1740/ 829-1772/ 891)  71.6 W48-W52 S15-S11
      9  319  599 (S17E59 L= 60)   770     97        2 ( 281/ 560- 354/ 639) -33.7 E61-E55 S19-S13
     10  771  382 (S35E27 L= 92)  1024    202        1 ( 718/ 316- 833/ 476) -85.6 E29-E24 S40-S27
     11 1199  453 (S34W06 L=125)   411    121        1 (1140/ 439-1260/ 468) -14.0 W01-W10 S35-S32
     12 1148  209 (S57E00 L=119)   768    104        1 (1098/ 191-1209/ 226)   4.8 E06-W06 S58-S54
  

    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