Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  923 1649 (N47W13 L=161)   354    116        0 ( 889/1605- 966/1691)  21.8 W08-W18 N45-N49
      2 1007 1384 (N28W08 L=156)   574    191        1 ( 917/1318-1067/1438)   9.5 W01-W13 N26-N30
      3 1148 1239 (N16W12 L=160)   209    105        0 (1096/1223-1190/1264)  -9.5 W09-W15 N15-N17
      4 1858 1104 (S12W59 L=207)  3927    482       10 (1785/ 873-1931/1261)  39.3 W47-W72 S22-S04
      5  702 1196 (N24E16 L=132)   871    219        1 ( 610/1137- 787/1252)   4.4 E22-E09 N22-N26
      6  416 1107 (N25E38 L=110)   464    108        0 ( 365/1064- 461/1144)  12.5 E42-E33 N23-N26
      7  516 1003 (N17E32 L=116)   272     88        0 ( 484/ 961- 546/1032) -68.2 E33-E30 N13-N18
      8 1560  993 (S08W31 L=179)   250     58        2 (1531/ 965-1588/1016) -45.0 W29-W32 S09-S06
      9 1486  814 (S16W22 L=170)  2465    529        5 (1299/ 702-1630/ 954)   8.7 W08-W34 S20-S11
     10 1438  872 (S12W20 L=168)   408     88        2 (1417/ 830-1456/ 915) -90.0 W18-W21 S14-S08
     11  525  601 (S07E41 L=107)   959    256        1 ( 418/ 532- 614/ 686) -50.7 E46-E36 S12-N00
     12 1543  611 (S31W23 L=171)   686    195        0 (1463/ 561-1611/ 679)  13.0 W16-W29 S32-S28
     13 1258  523 (S28W02 L=150)   290    105        0 (1244/ 466-1275/ 572)  63.4 W00-W03 S32-S25
     14 1146  447 (S30E08 L=140)  2199    501        6 (1054/ 314-1229/ 559)  41.2 E17-E01 S39-S23
     15  811  424 (S23E29 L=119)   622    164        2 ( 764/ 369- 850/ 497)  32.5 E35-E24 S26-S19
     16 1102  415 (S31E12 L=136)   513    171        0 (1024/ 391-1150/ 467)   0.0 E16-E06 S33-S28
     17  647  373 (S23E44 L=104)  1155    238        2 ( 591/ 296- 724/ 460)  85.9 E47-E40 S30-S16
     18 1414  364 (S46W09 L=157)   551     91        1 (1372/ 336-1454/ 393)   0.0 W05-W13 S48-S43
     19 1040  145 (S54E36 L=112)   979     91        2 ( 999/ 115-1079/ 175) -14.0 E42-E29 S58-S49
  

    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