Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  652 1680 (N40E42 L=124)   634    102        1 ( 604/1649- 686/1722) -15.9 E50-E36 N37-N42
      2  432 1617 (N32E60 L=105)   610     94        0 ( 416/1574- 452/1664) -28.6 E66-E54 N29-N35
      3 1279 1584 (N41W13 L=179)   755    151        1 (1224/1545-1344/1625) -29.1 W08-W18 N38-N43
      4  652 1436 (N23E30 L=135)   611    162        3 ( 618/1375- 701/1515)  83.7 E31-E28 N18-N28
      5 1503 1216 (N18W30 L=195)   228    123        0 (1480/1169-1538/1271) -56.3 W28-W32 N15-N21
      6  301  891 (S15E50 L=115)  4113    568        5 ( 190/ 658- 420/1084)  40.6 E75-E40 S31-S01
      7 1422  912 (S02W26 L=192)   245    139        0 (1364/ 871-1463/ 931) -15.9 W22-W29 S03-N00
      8 1702  742 (S10W51 L=217)   689    134        0 (1685/ 679-1731/ 805)  90.0 W49-W53 S13-S05
      9 1660  720 (S11W48 L=214)   335    115        0 (1625/ 681-1677/ 772)  80.5 W45-W49 S14-S08
     10 1252  746 (S14W17 L=183)   267     64        0 (1217/ 731-1293/ 762)   0.0 W14-W20 S14-S12
     11 1221  558 (S26W19 L=184)   396    124        1 (1164/ 529-1273/ 581)   7.1 W14-W22 S28-S24
     12 1641  509 (S25W55 L=221)  1155    113        1 (1588/ 480-1695/ 545)  24.0 W50-W61 S27-S22
     13 1270  435 (S34W26 L=192)   673    105        1 (1218/ 417-1328/ 456)  23.2 W22-W30 S35-S32
  

    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