Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1732 1516 (N32W60 L=107)   856    139        0 (1702/1485-1765/1539)  14.0 W55-W63 N30-N33
      2  377 1487 (N33E51 L=356)  1046    247        1 ( 329/1447- 438/1534) -19.7 E58-E44 N30-N36
      3 1547 1507 (N33W42 L= 89)   334     53        2 (1521/1487-1574/1534) -33.7 W39-W44 N31-N35
      4  702 1385 (N27E22 L= 25)   269     87        0 ( 659/1371- 749/1399)   8.1 E25-E18 N25-N27
      5 1112 1353 (N24W06 L= 54)   298     93        1 (1073/1324-1164/1376) -26.6 W03-W09 N22-N25
      6 1748 1263 (N16W54 L=101)   320     75        2 (1725/1217-1769/1298) -59.0 W51-W54 N13-N18
      7 1606 1245 (N16W40 L= 88)   384     89        0 (1572/1189-1632/1290) -71.6 W38-W42 N12-N18
      8  244 1222 (N15E58 L=349)  1073    184        2 ( 219/1134- 268/1288) -56.3 E62-E54 N10-N19
      9  891 1175 (N13E08 L= 39)   305     84        1 ( 847/1142- 932/1211) -33.7 E11-E05 N11-N15
     10 1286  969 (S00W16 L= 63)   249     88        0 (1241/ 950-1334/ 988)   0.0 W13-W19 S01-N01
     11 1182  569 (S26W10 L= 57)   282     85        0 (1137/ 559-1228/ 582)  -8.1 W06-W13 S26-S24
     12  401  552 (S27E48 L=359)   723    164        1 ( 364/ 523- 450/ 577)  21.8 E52-E42 S29-S24
  

    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