Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  862 1652 (N32E19 L= 94)   413     51        1 ( 831/1635- 901/1677)   0.0 E21-E15 N30-N34
      2 1372 1508 (N28W18 L=131)  1163    447        0 (1228/1436-1550/1629) -19.2 W07-W30 N25-N35
      3 1076 1440 (N20E01 L=112)   409    141        1 (1059/1369-1102/1499) -76.0 E02-W01 N15-N23
      4  924 1346 (N12E09 L=104)   483    179        0 ( 868/1280- 981/1415)  60.9 E11-E06 N08-N17
      5 1804 1168 (N13W52 L=165)  1590    271        4 (1757/1075-1877/1258) -55.0 W47-W60 N08-N18
      6 1864 1181 (N15W59 L=172)  1317    149        5 (1831/1131-1899/1246) -50.2 W54-W63 N11-N18
      7 1713 1110 (N08W44 L=157)  2081    407        3 (1577/ 943-1798/1245) -34.4 W32-W52 N00-N14
      8  254  879 (S20E54 L= 60)  4294    605        5 ( 128/ 754- 358/1060)  53.6 E69-E43 S27-S08
      9 1268 1018 (S03W14 L=128)   533    182        0 (1180/ 959-1320/1051) -12.5 W08-W17 S06-S01
     10  192 1026 (S11E59 L= 54)  1107    137        1 ( 120/ 987- 236/1050) -15.9 E69-E54 S13-S09
     11  540  874 (S19E29 L= 84)  1137    332        1 ( 466/ 811- 618/ 959) -53.1 E34-E24 S22-S13
     12 1838  731 (S12W64 L=177)  8499      0       14 (1621/ 539-1944/ 898)  56.8 W49-W77 S26-N00
     13 1219  839 (S14W14 L=127)   242     55        0 (1190/ 823-1255/ 858) -11.3 W11-W16 S15-S13
     14 1505  657 (S21W36 L=150)   788    190        1 (1473/ 586-1547/ 725) -37.9 W32-W41 S25-S16
     15 1672  678 (S18W49 L=162)   943    264        1 (1573/ 634-1770/ 721)  14.0 W40-W57 S21-S13
     16 1214  679 (S24W16 L=129)   230     79        0 (1181/ 647-1245/ 702)  45.0 W14-W17 S26-S22
     17 1407  630 (S24W30 L=143)   318     96        0 (1377/ 587-1428/ 676)  90.0 W28-W31 S27-S21
     18  451  548 (S40E41 L= 73)  1507    212        2 ( 359/ 518- 557/ 587)   7.8 E52-E30 S42-S37
  

    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