Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1563 1720 (N57W55 L=162)   791     90        1 (1528/1691-1594/1741)  13.0 W44-W70 N53-N60
      2  999 1393 (N24E06 L=102)   838    382        0 ( 920/1278-1121/1502) -40.9 E13-W03 N17-N30
      3 1503 1427 (N32W30 L=138)   499    130        1 (1449/1376-1537/1470) -50.2 W26-W33 N28-N34
      4  905 1312 (N17E11 L= 97)   510    163        2 ( 877/1227- 928/1377) -84.3 E12-E09 N11-N21
      5 1121 1312 (N20W03 L=111)   666    170        2 (1055/1249-1194/1344) -18.4 E01-W08 N16-N21
      6  998 1116 (N06E03 L=105)   190     63        1 ( 966/1101-1035/1130)  14.0 E04-W00 N05-N07
      7  824  763 (S17E10 L= 98)   491    112        2 ( 776/ 732- 885/ 792)  39.8 E12-E06 S19-S14
      8 1379  610 (S21W29 L=137)   401    123        1 (1353/ 573-1417/ 648) -18.4 W26-W32 S23-S18
      9  256  600 (S36E65 L= 43)   734     96        0 ( 221/ 567- 302/ 620)  12.3 E79-E56 S39-S33
     10 1435  390 (S36W43 L=150)   602    116        1 (1388/ 356-1487/ 428)  50.2 W39-W45 S38-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