Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1435 1638 (N45W34 L=263)   525     87        0 (1391/1612-1479/1656)  20.0 W28-W39 N42-N46
      2 1677 1603 (N42W62 L=291)  1621    113        3 (1630/1571-1740/1629)   2.1 W53-W81 N39-N44
      3  732 1544 (N35E25 L=204)   888    221        1 ( 676/1469- 789/1613) -37.6 E31-E18 N29-N39
      4  330 1426 (N24E57 L=172)   353     17        0 ( 310/1411- 354/1444)  -9.5 E59-E53 N23-N24
      5  698 1357 (N22E23 L=205)   267     57        1 ( 683/1322- 721/1387)  76.0 E24-E21 N19-N23
      6 1747 1131 (N11W51 L=279)  3508    397        8 (1656/ 959-1811/1294) -79.2 W43-W57 N00-N21
      7 1022 1199 (N13E01 L=228)   729    262        1 ( 934/1111-1094/1286) -48.0 E06-W03 N08-N18
      8  869 1205 (N13E10 L=218)   336     79        1 ( 850/1171- 892/1237) -63.4 E11-E08 N11-N15
      9 1256 1100 (N08W14 L=243)   394    150        0 (1230/1065-1318/1154) -33.7 W12-W18 N06-N11
     10 1086  981 (N00W04 L=233)   802    194        1 (1051/ 887-1127/1061)  84.3 W01-W06 S05-N05
     11 1363  825 (S08W22 L=251)  1091    252        1 (1326/ 716-1409/ 938)  81.9 W20-W24 S15-S01
     12  681  796 (S13E21 L=208)   255     74        0 ( 651/ 765- 710/ 829) -53.1 E22-E19 S14-S10
     13  945  808 (S11E04 L=224)   203     43        0 ( 911/ 793- 968/ 820) -14.0 E06-E02 S11-S10
     14 1869  727 (S14W71 L=299)  1441     81        4 (1845/ 685-1891/ 769)  90.0 W65-W77 S17-S11
     15 1763  728 (S14W56 L=284)   558     58        1 (1744/ 697-1792/ 751)  45.0 W53-W59 S15-S12
     16 1266  601 (S23W18 L=247)   641    122        0 (1217/ 558-1328/ 650)  45.0 W14-W21 S26-S19
     17 1214  493 (S31W16 L=245)   363     61        0 (1185/ 471-1241/ 523)  76.0 W13-W17 S32-S28
     18  937  230 (S55E04 L=224)  2346    265        4 ( 829/ 195-1073/ 260)   9.1 E17-W09 S60-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