Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1300 1810 (N51W54 L=171)  1569    207        2 (1208/1765-1393/1851)  -3.9 W39-W69 N48-N52
      2 1544 1701 (N36W64 L=181)   437     57        0 (1522/1677-1573/1717)   0.0 W60-W68 N34-N36
      3 1567 1362 (N16W43 L=160)   373     82        1 (1526/1339-1605/1387) -45.0 W40-W45 N13-N18
      4  334 1321 (N35E49 L= 68)   432    111        0 ( 305/1272- 363/1366)  56.3 E51-E46 N31-N37
      5 1690 1310 (N10W51 L=168)   533    133        1 (1637/1289-1756/1343)   0.0 W46-W58 N08-N11
      6  941 1240 (N21E01 L=116)   227    119        0 ( 916/1187- 979/1283)  51.3 E03-W01 N17-N22
      7 1041 1226 (N18W05 L=122)   400    201        0 ( 955/1198-1128/1251) -11.3 E01-W09 N16-N19
      8  676 1087 (N16E21 L= 96)   563    201        1 ( 627/1012- 741/1165) -65.6 E23-E17 N10-N21
      9 1803  973 (S12W53 L=170)  3930    432       13 (1736/ 794-1885/1106)  38.3 W46-W65 S21-S05
     10 1740 1048 (S06W48 L=165)   336     44        1 (1721/1016-1756/1073) -90.0 W46-W49 S07-S04
     11 1204  723 (S15W06 L=123)   758    254        3 (1172/ 592-1236/ 809)  69.0 W02-W08 S22-S09
     12  867  659 (S13E16 L=101)  1312    386        5 ( 769/ 575- 943/ 784) -65.0 E19-E12 S18-S03
     13 1117  739 (S12W01 L=118)   211     68        0 (1091/ 723-1157/ 761) -33.7 W00-W03 S13-S10
     14  497  700 (S05E40 L= 78)   276     90        1 ( 469/ 662- 533/ 729)  33.7 E42-E36 S07-S03
     15  390  545 (S14E54 L= 63)  1517    405        0 ( 335/ 400- 482/ 701)  87.4 E59-E49 S25-S03
     16  950  522 (S23E14 L=103)   408     91        0 ( 928/ 476- 972/ 571)  56.3 E16-E11 S26-S19
     17 1069  432 (S32E08 L=109)   725    118        1 (1012/ 405-1124/ 466)   7.1 E12-E04 S33-S29
     18  654  405 (S27E39 L= 78)   861    101        2 ( 619/ 348- 692/ 443)  76.0 E42-E35 S31-S23
     19 1331  280 (S52W09 L=126)   507     62        1 (1311/ 249-1358/ 305)  51.3 W07-W12 S55-S48
  

    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